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ABSTRACT 

Arcsecond-resolution X-ray imaging of the nucleus of the nearby starburst galaxy NGC 253 with 
Chandra reveals a well-collimated, strongly limb-brightened, kiloparsec-scale conical outflow from the 
central starburst region. The outflow is very similar in morphology to the known Ha outflow cone, 
on scales down to < 20 pc. This provides, for the first time, robust evidence that both X-ray and 
Ho; emission come from low volume filling factor regions of interaction between the fast energetic wind 
of SN-ejecta and the denser ambient interstellar medium (ISM), and not from the wind fluid itself. We 
provide estimates of the (observationally and theoretically important) filling factor of the X-ray emitting 
gas, of between ~ 4 and 40 per cent, consistent with an upper limit of ^ 40 per cent based directly on the 
observed limb-brightened morphology of the outflow. Only < 20 per cent of the observed X-ray emission 
can come from the volume-filling, metal-enriched, wind fluid itself. Spatially-resolved spectroscopy of the 
soft diffuse thermal X-ray emission reveals that the predominant source of spectral variation along the 
outflow cones is due to strong variation in the absorption, on scales of ~ 60 pc, there being little change 
in the characteristic temperature of the emission. We show that these observations are easily explained 
by, and fully consistent with, the standard model of a superwind driven by a starburst of NGC 253's 
observed power. If these results are typical of all starburst-driven winds, then we do not directly see all 
the energy and gas (in particular the hot metal-enriched gas) transported out of galaxies by superwinds, 
even in X-ray emission. 

Subject headings: ISM: jets and outflows 
X-rays: galaxies 



galaxies: individual (NGC 253) — galaxies: starburst 



1. INTRODUCTION 



Edge-on local starburst galaxies show unambiguous evi- 
dence for ^ 10 kiloparsec-scale, weakly-collimated, bipolar 
outflows (Heckman, Lehnert & Armus 1993). Current the- 
ory holds that these galactic superwinds are powered by 
the collective mechanical power of large numbers of Type 
II supernovae (SNe) and stellar winds, that result from the 
large population of massive stars formed in the starburst 
(Chevaher & Clegg (1985), henceforth CC). If this me- 
chanical energy is efficiently thermalized in the starburst 
region {i.e. converted back into the thermal energy of a 
hot gas by shocks, as SNe and stellar winds collide), then 
a pressure-driven outflow from the galaxy results. The 
hot gas blows out of the host galaxy's ISM along the mi- 
nor axis, forming the outflows seen in non-thermal radio 
emission, optical emission lines, and soft thermal X-ray 
emission in the halos of local starbursts. 

Superwinds are of cosmological interest as they trans- 
port large amounts of gas, in particular newly synthesized 
heavy elements, and energy, into the inter-galactic medium 
(IGM). Quantifying this mass, metal and energy transport 
in local starburst galaxies is essential for understanding 
the signiflcance of outflows from star-forming galaxies in- 
tegrated over the history of the Universe. However, even 



the basic physical properties of local superwinds such as 
mass outflow rates, energy content, abundances and kine- 
matics are uncertain. 

Measuring the physical properties of the hot gas driving 
these outflows is of crucial importance for several simple 
reasons. Firstly, the hot gas efficiently transports the ma- 
jority of the energy of the outflow (cf. Strickland & Stevens 
(2000)). Secondly, the SN-heated gas is thought to con- 
tain the majority of the newly synthesized metals. Finally, 
this energetic gas ultimately controls the ejection of mass 
from the galaxy (although the majority of the mass of the 
outflow may be in ambient ISM swept-up by the wind, this 
gas is accelerated to high velocity by the hot, fast, wind 
fluid, cf. CC). 

In principle. X-ray observations of the thermal emission 
from hot gas in superwinds can be used to measure the 
properties of the hot gas. Large amounts of ROSAT and 
ASCA X-ray data on starbursts already exists (see Read, 
Ponman & Strickland (1997); Ptak et al. (1997); Dahlem, 
Weaver & Heckman (1998), among others). However, the 
origin of the soft X-ray emission seen in local superwinds 
is still ambiguous. Plausible alternatives are: 

1. A hot, adiabatically-expanding, volume flUing fluid 
of stellar wind and SN-ejecta material, i.e. the 
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(a) 0.3-2.0 keV 




(b) 2.0-8.0 keV 



Fig. 1.— (a) Chandra ACIS S3 soft X-ray (0.3 - 2.0 kcV) image of the center of NGC 253. Tick marks correspond to 30". The image, 
shown on a square-root intensity scale, has been adaptively smoothed and background subtracted. Contour levels shown in black start at 
2.9 X 10^^ countss"^ arcsec"'^, and increase in steps of 3.2 X 10~^ countss"^ arcsec"^. White contours highlight the complex high surface 
brightness emission within the starburst region itself, at levels of (4.8, 6.5, 13.0, 26.0) xlO"* counts s"! arcsec . (b) The hard X-ray image 
(2.0 - 8.0 keV energy band) of the same region as in (a), again adaptively smoothed and shown using a square-root intensity scale. Contour 
levels have been chosen to highlight the clear lack of hard X-ray emission along the outflow lobes seen in (a). The black contour levels start at 
0.7 X 10"^ counts s-i arcsec (a factor 4 lower than the minimum used in (a)), and increase in multiples of two in surface brightness up to 
the minimum surface brightness level used in the contours of (a). In both panels the dashed and dotted lines show the position of the major 
and minor axes of the galaxy respectively, which intersect at the position of the brightest radio source (J2000.0, a = 00 47 33.14, S = -25 17 
17.2), which is presumed to be the nucleus (Ulvestad & Antonucci 1997). 



metal-enriched, energy containing "wind" itself (e.g. 
Fabbiano (1988); Bregman, Schulman & Tomisaka 
(1995)). 

2. A strongly "mass-loaded" wind, where the gas den- 
sity of the hot wind has been significantly in- 
creased by efficiently mixing in cooler ambient gas 
(e.g. Suchkov et al. (1996)), increasing the den- 
sity throughout most of the wind's volume. Mass- 
loading probably occurs primarily within the star- 
burst region. 

3. Emission from low filling factor regions of interac- 
tion between the hot, high velocity, wind and clumps 
or clouds of ambient ISM overrun by the wind. The 
X-ray emitting gas in and/or around the clumps 
may be heated by shocks (cf. CC), by localized hy- 
drodynamical mixing of hot wind and cooler cloud 
material, or by cloud evaporation due to thermal 
conduction. 

In the first two options X-ray observations probe the 

wind itself, i.e. the gas filling the majority of the volume, 
containing the majority of the energy, and the majority of 
the metals (although heavily diluted in case 2 by ambient 
ISM). In the third case X-ray observations probe low filling 
factor material that contains little energy, mass or metal- 
enriched gas. It is therefore essential to know the origin of 
the X-ray emission, both to correctly derive the properties 

^We assume the distance of 13 = 2.6 ± 0.7 Mpc to NGC 253, given 



of the X-ray emitting gas, and by extension those of the 
superwind as a whole. 

In this Paper, we demonstrate the first direct obser- 
vational evidence that soft thermal X-ray emission from 
starburst-driven superwinds arises in regions of interac- 
tion between denser ambient ISM and the hot, fast, wind 
fluid (case 3 above), based on Chandra ACIS X-ray obser- 
vations of the nearby starburst galaxy NGC 253. 

NGC 253 is one of the two brightest, closest^ and best 
studied starburst galaxies (the other is M82). The star- 
burst in NGC 253 has a bolometric luminosity of about 
2 X lO^°Z/0, and so is situated right at the "knee" of the 
starburst luminosity function (Soifer et al. 1987). Thus, 
starbursts with luminosities similar to NGC 253 dominate 
the starburst "emissivity" of the local universe. NGC 253 
is seen almost edge-on, ideally oriented for studying the 
superwind as it flows out into the galaxy halo. 

While both of the above attributes are shared by M82, 
NGC 253 offers two distinct advantages. Firstly, it lies be- 
hind a much smaller Galactic Hi column than M82 (10^° 
vs. 4 X 10^*^ cm~^ respectively). Since the superwind phe- 
nomenon is most conspicuous in soft X-rays, which are 
easily absorbed, this is a significant advantage. Secondly, 
NGC 253 (a late-type X* spiral) is much more typical of 
luminous starbursts than is M82 (a peculiar dwarf amor- 
phous galaxy). Thus, we expect the lessons learnt from 
NGC 253 to be more reliably generalizable. 



by Puche & Carignan (1988), for all calculations. 
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Fig. 2. — Logarithmically scaled continuum-subtracted Ha im- 
age of the same region of NGC 253 as shown in Fig. 1. The 0.3-2.0 
keV X-ray emission is shown as white contours. There is a clear 
similarity between the soft X-ray and Ha emission, in particular 
the southern outflow cone is strongly limb-brightened in both soft 
X-ray and Ha emission. 

2. DATA ANALYSIS 

Guest observer observations of NGC 253, using the 
AXAF CCD Imaging Spectrometer (ACIS^) on board the 
Chandra X-ray Observatory, were obtained on 1999 De- 
cember 16, with the nuclear star burst region placed on 
the back- illuminated CCD chip S3. 

Here we report on the soft X-ray properties of the out- 
flow within the central kiloparsec-scale region of NGC 253, 
and defer further discussion, and detailed spectroscopic 
study, of X-ray emission from point sources and the larger- 
scale diffuse X-ray emission to forthcoming papers (Weaver 
et al. (2000); Strickland et al. (2000)). 

Data reduction and analysis was performed using CiAO 
(vl.1.4) and Heasoft (v5.0.1). The data were repro- 
cessed using the latest gain correction file available at the 
time of writing (acisD1999-09-16gainN0003.fits), screened 
to remove "bad" events, and time filtered to remove pe- 
riods of high background ("flares," where the total count 
rate in any chip deviated by more than 5(t from the mean) . 
This resulted in a remaining exposure length of 12862s on 
chip S3. 

Background subtraction for both imaging and spec- 
tral analysis was performed using the background event 

files and software provided by the Chandra X-ray Center 
(CXC), scaled to the exposure time of our observation. 

Spectra were extracted using Pulse Invariant (PI) data 
values, in order to account for the spatial variations in gain 
between different chip nodes. In order to use as the fit 
statistic the spectra were rebinned to achieve a minimum 
of at least 10 counts per bin. Given the calibration un- 
certainties in the ACIS BI chip spectral response (in par- 
ticular at energies below E ^ 1 keV), we seek only to use 



spectral fitting to qualitatively characterize the emission, 
and do not place too much signiflcance in the numerolog- 
ical results of the fits themselves. 

Spectral responses were created using the latest calibra- 
tion files available for observations taken with the CCD 
temperature of -110 °C (the April 2000 release), following 
the example procedures provided by the CXC"*. We used 
the spectral response appropriate for the fiux-weighted 
center of the spatial region from which any spectrum was 
extracted (typically a region ^ 40 x 100 raw pixels in size, 
see § 4). This does not fully account for the spatial vari- 
ation in the ACIS spectral response, but is a good first 
approximation, given the current lack of publically avail- 
able software to construct spectral responses for spectra 
from regions largc;r than the 32 x 32 pixel regions over 
which a single response strictly applies. 

The data products initially released to us by the CXC 
suffered from processing errors that introduced an ^ 9" 
shift in the absolute coordinate system for the data. For 
the purposes of this Paper, we registered the Chandra data 
against ground-based R-band and continuum-subtracted 
Ha data, based on 9 secure optical cross-identifications of 
X-ray sources seen in the Chandra observation. This is 
accurate to ~ 0.9" (root mean square). 

3. GENERAL FEATURES OF THE X-RAY DATA 
3.1. Soft X-ray emission 

Many different sources of X-ray emission are visible in 
the Chandra soft X-ray (0.3 2.0 keV energy band) image 
of NGC 253 (Fig. la). Emerging from the nuclear star- 
burst region, bright diffuse soft X-ray emission can be seen 
extending away from the plane of the galaxy to the south 
east and north west for a distance of ^ 1'.2 and 0^6 respec- 
tively (~ 900pc and ~ 450pc, assuming D = 2.6 Mpc). 
This is the innermost part of the starburst-driven super- 
wind, which has been traced with RO SAT out to ^ 9kpc 
from the nucleus (Dahlem, Weaver & Heckman 1998). 
The southern lobe seen by Chandra is well described as a 
limb-brightened conical outflow, of opening angle 9 ~ 26°, 
with a truncated base of diameter ~ 24", similar in size 
to that of the starburst region itself The outflow cone 
is also seen in Ha images, and the positions and size of 
the X-ray and Ha outflows match up almost exactly (see 
Fig. 2). Beyond a distance of ~ 600 pc from the nucleus, 
the X-ray emission becomes too faint to trace the limb- 
brightening further, and the Ha morphology at this point 
becomes confused due to emission from what are known 
to be bright Hii regions along the spiral arm seen in the 
background (McCarthy, Heckman & van Bruegel 1987). 
We analyze the spatial structure of the limb-brightened 
region of the southern outflow lobe in more detail in § 5. 

Due to NGC 253's inclination, the northern outflow lobe 
is seen through the disk of the galaxy. Strong absorption 
(see § 4) makes the inner region of the northern outflow 
less obviously conical, and it can not be traced as far as 
the southern outflow cone in the Chandra data. 

Many X-ray point sources are also clearly visible, pre- 
sumably high mass X-ray binaries or supernova remnants. 

The peaks in the X-ray surface brightness at the south- 



'^Chandra Observatory Guide http;//asc. harvard.edu/udocs/docs/docs. html, section "Observatory Guide" , "ACIS" 
^See the ClAO Science Threads: http://asc.harvard.edu/ciao/threads/threads.html 
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eastern-most extent of the limb-brightcned region are re- 
solved structures, and are not point-like. The effective 
spatial resolution of the data within the region shown in 
Fig. 1, measured by fitting 2-dimensional Gaussian mod- 
els to the brightest point sources, is ^ 0'.'9. As the pixel 
size of the ACIS instrument is 0'.'492, point sources are 
easily distinguished from diffuse emission in unsmoothed, 
unbinned, images, as most of the counts of a true point 
source can be found within a few adjacent pixels. 

3.2. Hard X-ray emission 

Hard X-ray emission from NGC 253, in the energy band 
2.0 - 8.0 keV, is dominated by point source emission (see 
Fig. lb). Several point-like sources are visible in the cen- 
tral starburst region (at least three within the central 8" 
100 pc] radius), although these are not the brightest 
X-ray sources in NGC 253. 

The nucleus is surrounded by what appears to be ex- 
tended hard X-ray emission, roughly spherical but some- 
what elongated preferentially along the major axis of the 
galaxy, which wc can trace out to a radius of ^ 8" (100 
pc) with confidence (in that this size is apparent to the eye 
in the raw, unsmoothed, data). At lower surface bright- 
ness levels the increasing apparent size and sphericity of 
this diffuse hard X-ray emission may be an artefact of the 
adaptive smoothing algorithm used. The origin and prop- 
erties of this emission will be discussed in Weaver et al. 
(2000). 

There is no significant hard X-ray emission from the the 
southern outflow cone, which is so clearly seen in the soft 
X-ray band. We shall discuss what this implies for the 
presence of very hot gas {kT > 2 keV) in the outflow in 
§ 4 & 5. 

4. SPECTRAL STRUCTURE WITHIN THE OUTFLOW 

We shall not attempt a rigorous, quantitative, spectro- 
scopic study of the outflow in this Paper given the uncer- 
tainties in the spectral calibration of the ACIS instrument. 
We believe that the data can support two robust conclu- 
sions: 

1. The primary cause of spectral variation along, and 
between, the northern and southern outflow cones 
is due to variations in the absorption column, and 
not due to significant temperature variations along 
the outflow. 

2. Features in the ACIS spectra of the diffuse emis- 
sion show clear evidence for X-ray lino emission 
from highly ionized a-elements (e.g. Mg and Si, see 
Fig. 3). This clearly demonstrates the thermal na- 
ture of the X-ray emission. 

4.1. Simple spectral fitting 

To roughly characterize the diffuse X-ray emission from 
the inner kiloparsec of the superwind, i.e. an approximate 
temperature and amount of foreground absorption, we di- 
vided the outflow lobes and nuclear starburst region into 
the set of regions shown in Fig. 3. From each region we 
then extracted the spectrum of the difi^use emission, hav- 
ing excluded any detected point sources, and used a simple 
one-temperature hot plasma model to fit the spectrum. 



The region of the southern outflow lobe most clearly 
limb brightened, in both the X-ray and Ha data, falls 
within our spectral regions SI and S2. As this limb- 
brightened region if of particular interest (and is the focus 
of § 5), we also extracted and fitted the spectrum (Fig. 4) 
of this region (henceforth called the CLE region). 

Only the data within the energy range between E = O.'d 
to 3.0 keV was used in the spectral fitting, given that the 
2.0 - 8.0 keV ACIS count rates for all the regions (ex- 
cept the starburst region NC and the northern lobe region 
Nl) are essentially zero (see Table 1). The lack of hard 
X-ray emission in the southern outflow lobe is graphically 
demonstrated in Figs. 1 & 4. 

To begin with, we fitted a variable abundance, single 
temperature Mekal hot plasma model to the combined 
spectra from the southern outfiow cone, (regions SI, S2 & 
S3 in Fig. 3) where a hardness ratio analysis (see § 4.2) 
already indicated low absorbing columns and similar spec- 
tral shape. The abundance ratios (with respect to the So- 
lar value) of the a-elements were constrained to be equal 
to each other, while the abundance ratios of the other el- 
ements were forced to be equal to that of Iron. 

TABLE 1 

Spectral pits to dippuse X-ray emission regions 



Region Count rate" N^'' 


kT'' 




0.3-2.0 keV 2.0-8.0 keV 10^° cm-^ 


keV 





Nl 

NC 

SI 

S2 

S3 

CLE 



0.64 ±0.16 


92.2t\ll 


0.59; 


1.08 ± 0.17 




0.66: 


0.15 ±0.16 




0.63: 


0.07 ±0.16 


< 3.4 


o.5i: 


0.05 ± 0.16 


< 5.0 


0.46: 


0.09 ± 0.16 




0.57: 



+0.10 
-0.08 
+0.05 



+0.13 



128.5 / 67 



" Background-subtracted ACIS S3 count rates in units of 
10~^ count s~^. 

* Confidence regions are 90% confidence for three free parameters 
(Ax^ = 6.25). 



TABLE 2 

Hardness ratios in region of clear limb brightening 



Region 






Count 


rates'" 












deg 


0.3 - 


1.0 keV 


1.1 - 


2.0 keV 








Entire CLE 


130- 


-160 


6.06 


±0.22 


0.87 


± 0.09 


-0.75 


± 


0.04 


East limb 


130- 


-137 


1.58 


±0.11 


0.21 


±0.07 


-0.77 


± 


0.09 


Center 


137- 


-144 


0.72 


±0.08 


0.10 


±0.06 


-0.76 


± 


0.16 


West limb 


144- 


-160 


3.65 


±0.17 


0.55 


±0.07 


-0.74 


± 


0.06 



" Range of position angle in the azimuthal region (see Fig. 6. 

Background-subtracted ACIS S3 count rates in units of 
10~^ count s~^. 

° Hardness ratio Q = {H — S)/(H + S), where H is the count rate 
in the 1.1 - 2.0 keV band, and S is the count rate in the 0.3 — 1.1 
keV band. 



This spectral fit, although of a region of very low 
and relatively uniform hydrogen column density (to mini- 
mize possible spectral complexity), gives strongly sub-solar 
abundance ratios {Zo, ~ 0.30tg;i^ Zq, Zpe ~ 0.131^;^^ Zq, 
90% confidence in one interesting parameter). This is 
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Northern outflow region N1 




chonnel energy (keV) 



Fig. 3. — Three color, logarithmically-scaled, soft X-ray image of the central 2' X 2' (~ 1.5 X 1.5 kpc^) region, showing the spectral variation 
along the outflow lobes. Emission between 0.3 < E{keV) < 0.6 is shown in red, emission between 0.6 < E(keY) < 1.1 in green, and emission 
between 1.1 < i?(keV) < 2.0 in blue. The south-eastern outflow cone, on the near side of the disk, becomes progressively "soft" as it rises 
above the absorbing Hi gas in NGC 253's disk and starburst region. In contrast, only the highest energy X-ray photons from the north-western 
lobe of the outflow make it through the disk without absorption. The white dashed boxes show regions from which spectra of the difi^use 
emission, shown in the panels on the right, were extracted. The yellow dashed line outlines the region most clearly limb brightened in both 
X-ray and Ha images. Dashed red circles show regions of emission from point sources that were excluded in extracting the diffuse spectra. 
The spectra of the nuclear region and the northern outflow clearly show signiflcant absorption at photon energies E < 1 keV. 



Hkely to be a combination of further spectral complexity 
on small scales (emission from multi-phase gas, cf. Strick- 
land & Stevens (2000) and Weaver, Heckman & Dahlem 
(2000)) and complications arising from the uncertain and 
spatially varying instrument response (e.g. we believe the 
failure of the model to fit the spectrum of region NC in 
the energy range 1.6 - 2.0 keV to be due to problems in 
the spectral response used), rather than true evidence of 
low metal abundances in the gas. 

This spectral model was used as a template fit to all 
the separate regions along the outflow lobes (see Fig. 3), 
allowing only absorption column density, gas tempera- 
ture and absolute model normalization to vary (see Ta- 
ble. 1). Best fit hydrogen columns vary by a factor ~ 45, 
from TVh « (9 ± 1) x 10^^ cm~^ (Northern region Nl) to 
Nn « 2 X 10^°cm~^ (Southern region CLB). In contrast 
characteristic temperatures vary by less than a factor of 
2, from kT « 0.66lJ] Jg keV (NC) to kT w 0.46l[5 iJ keV 
(S3). 

4.2. Hardness map 

To fully exploit Chandra^s spectral-imaging capabilities 
we constructed the "hardness map," shown in Fig. 3, of the 
central ~ 1.5 kiloparsec square region of NGC 253. This is 
a three color image of the data between (approximately) 
0.3 - 0.6 keV, 0.6 - 1.1 keV and 1.1 - 2.0 keV in energy. 



The southern outflow cone is clearly much softer (a higher 
flux of low energy X-ray photons) than the diffuse emis- 
sion from the starburst region, which itself is softer than 
the emission seen from the northern outflow on the far 
side of the disk, consistent with the simple spectral fitting 
discussed above. 

These three energy bands were chosen based on the 
observed spectral shape of the diffuse emission spectra. 
Emission in the 0.3 - 0.6 keV energy band, relatively 
prominent in the SI to S3 and CLB regions, is the first 
to be absorbed away as the foreground absorption column 
increases. The relative prominence of the 0.6 - 1.1 keV 
hump (which contains the unresolved Fe L shell complex 
of lines) in the nuclear region (region NC) is largely due to 
the higher column in this region having removed much of 
the emission below E = 0.6 keV. In region Nl the column 
is high enough to absorb away much of emission between 
0.6 - 1.1 keV, making the 1.1 - 2.0 keV energy band ap- 
pear much stronger (and hence the spectrum harder) than 
in the other spectra. As the spectral fits in Table 1 show, 
the variation in spectral hardness of the diffuse emission 
is mainly due to increasing absorption. 

Note that the nucleus of the galaxy lies near the north 
western edge of region NC (compare Figs 1 & 3), so that 
region Nl partially covers the most obscured part of the 
starburst region as well as the innermost part of the north- 
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ern outflow lobe. 

The measured in-flight spectral resolution of the ACIS- 
S3 chip is ^ 0.12 keV (full width at half maximum), at 
the operating temperature of -110 °C that applied to these 
observations (sec the Chandra Observatory Guide). The 
three energy bands used in Fig. 3 arc therefore largely in- 
dependent of each other (as can bo scon from the spectra), 
i.e. "cross-talk" between the bands is minimal. 

These Chandra observations show that emission from 
the wind is spectrally complex, due to spatial variations 
in absorption over scales as small as ~ 5" (~ 60 pc), possi- 
bly with multiple temperature components on even smaller 
scales. The southern outflow cone is relatively spectrally 
uniform, having risen above the location of absorbing gas 
in NGC 253'h disk. Without high spatial resolution. X- 
ray spectroscopy of diffuse X-ray emission in star-forming 
galaxies will continue to yield ambiguous results due to 
mixing regions of different temperature and/or absorbing 
column together, in addition to problems due to blending 
unrelated point sources into the spectra. 




0.5 1 2 

channel energy (keV) 



Fig. 4. — Chandra ACIS S3 chip X-ray spectrum of the Umb 
brightened region of the southen outflow lobe (region CLB) , fit- 
ted with a single temperature Mekal hot plasma model. Note the 
lack of any additional hard X-ray component at energies above 2 
keV. See § 4.1 and Table 1 for details of the spectral fit. 

5. SPATIAL STRUCTURE OF THE OUTFLOW 

The striking similarity between X-ray and Ha emission 
in the southern outflow cone can be used to constrain the 
origin of the soft X-ray emission seen in the superwind, and 
its physical properties. Although it had been known from 
even early Einstein X-ray observations that the general 
spatial extent and distribution of X-ray-emitting gas in 
NGC 253 and M82 was similar to the Ha emission (Wat- 
son, Stanger & GrifSths (1984); Fabbiano (1988)) the 
^ 5" spatial resolution, the 3 5" positional imcertain- 
ties, and limited sensitivity, of both Einstein & ROSAT 
High Resolution Imagers (HRIs) only implied that X-ray 
and Ha emission were similar over scales > 60 - 100 pc. 
In contrast, our Chandra observations can be used to in- 
vestigate the X-ray / Ha correlation down to the effective 
resolution of the data (~ 0'.'9 w 11 pc). 

5.1. Azimuthal surface brightness profiles 

We looked at both azimuthal surface brightness profiles 
of the limb-brightened region, and at surface brightness 



"slices" constructed perpendicular to the minor axis of the 
galaxy. Both methods give similar results, but azimuthal 
profiles have the advantage of preserving the conical, di- 
verging, structure of the outflow. We use the degree of 
limb-brightening to constrain the fraction of X-ray and 
Ha emission that arises in the shell or wall of the outflow, 
compared to emission from within the central core of the 
outflow itself. 

We created azimuthal profiles of background-subtracted 

X-ray and Ha surface brightness, from the region of clear 
limb-brightening (see Fig. 5), using 1.0° and 0.5° az- 
imuthal bins respectively, taking care to exclude regions 
of emission from unrelated point sources. The resulting 
surface brightness profiles are shown in Fig. 6. 

Both X-ray and Ha azimuthal profiles peak at the same 
position angle on the western limb of the outfiow. In both 
cases the eastern limb of the outflow cone is ~ 0.75 the 
intensity of the western limb, and the surface brightness 
decreases inwards of the eastern limb to a minimum value 
of ~ 0.3 - 0.4 of the peak brightness (although the X-ray 
peak is slightly offset to the interior of the cone from the 
Ha). There is not such a clear decrease in brightness in- 
wards of the western limb, but the images (Fig. 5) show 
additional "clouds" apparently within the outflow in both 
X-ray and Ha images. These features might even be unre- 
lated emission associated with the underlying disk, rather 
than part of the outflow itself. 

We use a simple model to quantify the fraction of emis- 
sion coming from within the core of the outflow, and what 
fraction comes from the apparent shell. We assume a con- 
ical outflow, with a uniform emissivity core /core of open- 
ing angle 0core, surrounded by a shell of emissivity /shell 
and outer opening angle ^^sholl (sec Fig. 5). Predicted az- 
imuthal proflles are convolved with Gaussian mask, to take 
account for the effective spatial resolution of the X-ray and 
Ha data (FWHM = 0'.'9 [= 0.6° in position angle in the 
azimuthal profiles] for the X-ray and 1'.'5 for the Ha data), 
and compared by eye with the data. We primarily attempt 
to match the position, width and intensity of the eastern 
limb, given the more confused state of the western side 
of the profiles. The relative intensity of the eastern limb 
with respect to the central surface brightness minimum 
then constrains the emissivity of any material within the 
core of the outflow. The best flt-by-eye shell models are 
shown in Fig. 6. Attempting to better match the inten- 
sity of the western limb would increase the peak-to-core 
intensity contrast and lead to lower values of /core//sheii- 

Both the X-ray and the Ha data are consistent with a 

hollow (/core //shell = 0) COUical OUtfloW with ^shell = 26° 

and ^corc ~ 20°. The maximum value of /coro//shoii con- 
sistent with the data is ^ 0.2 (X-ray) or ^ 0.1 (Ha). Both 
the X-ray and Ha emission must therefore arise predomi- 
nantly in a thin boundary layer, between a faster and/or 
hot tenuous wind of SN-ejecta (that invisibly fills the core 
of the outflow), and dense cool ambient ISM surrounding 
the walls of the cavity. 

We looked for spectral variation between the center and 
walls of the southern outflow cone, in an effort to detect 
hotter gas from within the cone, but foimd no statisti- 
cally significant differences. We divided the limb bright- 
ened region into three sectors, corresponding to the east- 
ern limb, the low surface brightness core, and the western 
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Fic;. 5. — (a) Spatial model assumed for the origin of the X-ray and Ha emission, (b) Chandra 0.3 — 2.0 kcV image, with the region of clear 
limb-brightening between i?min and i?max, used in extraeting the azimuthal surface brightness profile, shown by the curved dashed lines. This 
region lies between 27" and 47" (~ 340 - 590 jjc) from the major axis of NGC 253. These lines, along with the dashed lines parallel to the 
outflow, delimit the spatial region used to calculate gas masses, densities and filling factors as described in § 5. (c) As (b), but showing the 
same region on the Ha image. 



hmb, based on the azimuthal profile shown in Fig. 6. As 
there were too few counts in these sub-regions for spectral 
fitting wc constructed a hardness ratio Q from the back- 
ground subtracted count rates in the two energy bands 
0.3 - 1.1 keV and 1.1 - 2.0 keV, and calculated Q for 
the entire limb brightened region, as well as the three sub- 
regions (see Table. 2). Simulations with Xspec had shown 
this choice of energy bands to be reasonably sensitive; to 
changes in spectral hardness for thermal emission with 
temperature fcT ~ 0.1 to 3.0 keV and hydrogen columns 
of between a few x 10^° to a few x 10^^ cm^^ . This hard- 
ness ratio is more sensitive to temperature variation than 
using a simple ratio of 0.3 - 2.0 keV verses 2.0 - 8.0 keV 
count rates, for example. Both the limbs and core of the 
limb brightened region have hardness ratios statistically 
indistinguishable from one another or from the region as 
a whole. This is consistent with the majority of the emis- 
sion, even within the pro jected center of the outflow cone, 
being from the bright walls of the cavity. Nevertheless, 
the uncertainty in the hardness ratio of the core region is 
large enough to correspond to an uncertainty in tempera- 
ture (from the nominal kT ~ 0.6 keV) of approximately a 
factor two. 

The spectral uniformity of the southern outflow cone al- 
lows us to interpret variations in surface brightness directly 
as variations in emission integral (EI = JncnudV). At- 
tempting to explain the central surface-brightness decre- 
ment of 50 - 60%, with respect to the apparent walls 
of the cone, as being purely due to absorption would re- 
quire an order-of-magnitude increase in absorption column 
in front of the center of the cone. This would also make 
the surprising degree of similarity between the Ha and X- 
ray azimuthal profiles an astounding coincidence, which 
seems unlikely. For example, assuming emission from a 
hot plasma with parameters equal to the best-fit model for 
the CLB region (see Table. 1) for the "walls" of the cavity, 
then reducing the fiux in the center of the cone to 40% of 
its original value by extra absorption requires increasing 
the column by a factor 12 to Nu ^ 2.5 x 10^^ cm~^ . This 
would lead to a strong decrement in the ratio of 0.3 - 0.6 
keV to 0.6 - 1.1 keV count rates in the center of the cone. 



relative to the walls of the cone, which is not seen. Con- 
structing a hardness ratio from these two energy bands, 
in a manner analogous to those shown in Table 2, gives 
a hardness ratio of 0.30 ± 0.12 for the center of the cone, 
compared to a value of 0.35 ±0.04 for the entire CLB region 
and the predicted value of 0.73 for A^h = 2.5 x 10^^ cm~^ . 

The inferred opening angle of the outflow, ~ 26°, is 
very similar to the 9 ^ 30° opening angle for the out- 
flow in M82. The inferred diameter of base of the out- 
flow, i.e. in the starburst region, is 24" (300 pc) which 
agrees very well with the ^ 20" diameter region over which 
young SNRs are seen in radio observations (Ulvestad & 
Antonucci 1997), consistent with the idea of a SN-driven 
outflow. 

5.2. Derived physical properties of the outflow 

The maximum volume filling factor of the X-ray & 
Ha emitting material in the outflow cone is then r] ~ 0.4, 
based on the assumed geometry. If both X-ray & Ha- 
emitting gas arise in a shocked or conductively heated 
boundary layer, in close proximity to each other, then 
their pressures should be similar on local scales {i.e. < 
a few 100 pc. On larger scales the outflow is driven by 
the over-pressure of the starburst region compared to the 
ambient ISM or IGM). We can then infer the true density 
of the X-ray emitting gas and its filling factor. In a 30°- 
wide azimuthal region centered on the position angle of 
the outflow cone and between the radii defined in Fig. 5, 
the emission integral EI w 5.5 x 10^" cm"'' . This is based 
on the observed 0.3 - 2.0 keV Chandra ACIS-S count rate 
of 0.070±0.002 counts s"^ in this region, and assuming for 
simplicity emission from a fcT = 0.5 keV hot plasma of So- 
lar abundance, absorbed only by the Galactic foreground 
hydrogen column of A'h — 10^° cm"^ . The total volume 
{i.e. shell -|- core) of the conical frustum defined by our 
assumed geometry is F = 1.4 x 10^^ cm^, which then gives 
ric ~ 6.2 X 10~^?7"^/^ cm"'' . From McCarthy, Heckman 
& van Bruegel (1987), Pna/k - 1 to 4 x lO'^Kcm^^ in 
the Ha emitting gas, in the region between 250 and 500 pc 
from the nucleus. The true density and volume filling fac- 
tor of the X-ray-emitting gas is then rie ~ 0.1 to 0.3 cm~^ 



8 



Strickland et al. 



and 77 ~ 40 to 4 per cent (cquivalcntly ~ 100 to 10 per 
cent of the geometrical volume of the shell). 

The main uncertainty affecting the derived EI is the as- 
sumed metal abundance. As we believe the fitted metal 
abundance to be unrealistically low, for the reasons given 
in § 4.1, we prefer to assume a plausible model for the 
emission instead of relying too heavily on simplistic spec- 
tral fits. Assuming the gas temperature is kT ~ 0.6 keV 
instead of 0.5 keV only alters the EI by ~ 0.5 per cent. 




120 130 140 150 160 170 



Position angle (degnees) 

Fig. 6. — (a) Normalized azimutlial surface brightness profiles 
of the southern outflow cone, in Hq emission (diamonds, dashed 
line) and 0.3 — 2.0 keV X-ray emission (erosses, dotted line), (b) 
Hq profile, with best-fit hollow cone model (/core //shell = 0) over- 
laid (solid line). A model with /core//shell = 0.2 (dashed line) 
over-prediets the central surface brightness between position an- 
gles 134 - 142°. (c) As (b), but for the soft X-ray emission. 
Models with /core//shell > 0-3 (dashed line) fail to fit the data 
between position angles 134 — 142° . 

The total mass in X-ray emitting gas, within this spe- 
cific region, is Mx-ray ~ 1.5 - 6 x 10^ M©. We used the 
spectroscopic observations of Keel (1984) to flux cahbrate 
our Ha image (assuming a constant Ha/ [Nil] ratio), and 
hence calculate the Ha flux and gas mass in the same con- 
ical frustum as used in the X-ray analysis. Assuming the 



temperature of the Ha emitting gas is Thc = lO^'K, and 
that the electron number density of this gas is in the range 
He = 100 - 25 cm~^ (based on McCarthy, Heckman & van 
Bruegel (1987)), the mass of the Ha emitting gas in this 
region is very similar to the mass of X-ray emitting gas, 
~ 1 ~ 5 X 10'' Mq. X-ray and Ha gas masses are 
similar, despite the large difference in gas density, as the 
derived filling factor of the Ha emitting gas within this 
region is substantially smaller than the filling factor of the 
X-ray gas (?7Ha 0.4 7 x 10^'*). 

Low X-ray filling factors are a prediction of current 
hydrodynamical models of superwinds. In a parameter 
study of superwind models (Strickland & Stevens 2000), 
based roughly on M82's superwind, the globally-averaged 
filling factor of the X-ray emitting gas in the ROSAT 0.1 
- 2.4 keV energy band (roughly equivalent to the 0.3 - 
2.0 keV ACIS-S energy band we have used) was always 
< 2 per cent. 

Is it possible to have a wind powerful enough to com- 
press the ambient ISM to pressures P/k > lO^K cm~^ , 
yet also tenuous enough to be largely invisible in X-ray 
emission? In short, the answer is yes! 

The maximum emission integral of X-ray-emitting gas 
not from the shell, based on /coreZ-^shcii — 0.2, is Elmore = 
1.2 X 10^° cm-3 {kT = 0.5 keV) or 3.7 x 10^° cm-'"* {kT = 
2.0 keV, a realistic value for the wind fluid temperature 
in this region based on the CC model). These numbers 
correspond to - 0.015 ACIS S3 count s^^ in the 0.3 - 2.0 
energy band, and imply count rates in the hard 2.0 - 8.0 
keV energy band of ~ 3 x 10~^ (assuming kT = 0.5 keV) 
or 1.6 x 10-3 count s"^ {kT = 2.0 keV). The total ob- 
served 2.0 - 8.0 keV count rate in this region is only 
(0.9 ± 1.6) X lO^'' count s^^, which can be accounted for 
purely by the kT = 0.6 keV emission from the shell (ex- 
trapolation from the 0.3 - 2.0 keV count rate predicts a 
2.0 - 8.0 keV count rate of - 1.4 x 10"^ count s'^). This, 
and the hardness ratios given in Table 2, do not rule out a 
spectrally harder core component to the X-ray emission in 
the limb brightened region, given the poor statistics. Nev- 
ertheless they suggest that the upper limits on the core 
emission integrals given above are conservative. 

In the CC model of a freely-expanding wind, ram pres- 
sure dominates over thermal pressure at distances beyond 
a few starburst-region radii. For our model to work we 
require a wind ram pressure such that Pram = Pw^^ > 
-Px-iay ^ Piia, as it is the wind that shocks and com- 
presses the ambient ISM. This is a lower limit — the 
wind ram pressure Pram can be several times Pua as the 
Ha emitting material is also being swept-up and acceler- 
ated to velocities of several hundred kilometers per sec- 
ond by the wind. Taking this into account, let us assume 
Pram = V'-Phq- luspcction of our existing hydrodynamical 
simulations of superwinds shows 1 < V ^ 10) there being 
a large scatter in ^ from region to region within the wind. 
The lower limit on the emission integral of the wind fluid, 
required to give acceptable wind ram pressures in the same 
region as the limb-brightened Ha and X-ray emission, is 

-B/w > V{1 - ryx-ray " Vila) X {ij^Plia / ^fj/mnvl)'^ ~ 1 - 

13 x lO^^V'^^^a'^ cm""^ , where V3 is the wind velocity in 
units of 3000 km s~^ (^3 ~ 1 for standard mass and en- 
ergy injection rates, e.g. see Leitherer & Heckman (1995)), 
and assuming tp = 1 to give a lower limit. Assuming that 
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tp = 3 is typical of the wind on average, then the wind 
velocity must be V3 > 0.3 - 0.6, in order not to violate the 
observed upper limits on wind emission integral (using the 
lower, and more conservative, of the two values given in 
the previous paragraph), which seems entirely reasonable. 

These limits also make sense in terms of mass transport 
rates in the wind. The observed SN rate of ~ 0.05 yr^^ 
(Colina & Perez-Olea 1992) corresponds to an expected 
~ O.5M0yr~^, using the Leitherer & Heckman 
(1995) models, and assuming no additional mass-loading. 
Based on the pressure balance argument, the total mass 
loss rate (assuming two identical outflow cones) in the 
wind fluid alone is M„ > 0.032 - 0.13 i^v^'^ Mq yr^^ The 
observational upper limit on the emission integral of the 
core of the outflow gives Mw < 2.2?;3 yr~^, using the 
higher of the two EI core estimates. 

It is interesting to note that a clear detection of the wind 
component would have placed strong constraints on wind 
density, velocity, energetics and mass outflow rate. Sepa- 
rating the faint X-ray emission due to the wind from the 
brighter X-ray emission from wind/ISM interfaces in local 
starburst galaxies requires arcsecond spatial resolution — 
a job only Chandra can do. 

5.3. Additional discussion 

How robust are these conclusions on the origin and flll- 
ing factor of the X-ray emitting gas. given the limited size 
of the region of unambiguous limb-brightening, the sim- 
plicity of the spatial model assumed and the uncertain 
value of the gas emission integral? 

Although these factors probably make the estimated 
densities and filling factors of the X-ray emitting gas un- 
certain to an order of magnitude or so, the general pic- 
ture of the X-ray and Ha emission being dominated by 
moderate-to-low filling factor gas, arising in the interac- 
tion between the wind and denser ambient ISM, is likely to 
remain true. The simple analytical arguments presented 
in § 5.2 certainly demonstrate that our conclusions are 
physically plausible. 

Within the region we have considered, the filling fac- 
tor of the gas that is the dominant X-ray emitter must be 
less than unity, as it is impossible for a volume-filling gas 
to reproduce the observed limb-brightening without invok- 
ing fortuitous absorption along the center of the outflow 
cone. Such additional absorption would be detectable in 
the spectra and hardness maps, and is therefore inconsis- 
tent with the observed spectral uniformity of the southern 
outflow cone. 

That strong limb-brightening is only observed in a nar- 
row range of heights above the plane of the galaxy is not 
surprising. Within the immediate vicinity of the starburst 
region Ha and X-ray emission not directly associated with 
the superwind itself will hide any limb-brightened out- 
flow component. At much larger distances, as the surface 
brightness of the superwind decreases, limb-brightening 
would become increasingly difficult to detect. A well de- 
fined cavity or chimney structure with dense walls might 
only exist within a few scale heights of the disk. Further 
out in the wind, the dense gas in the outflow probably ex- 
ists as isolated dense clumps and clouds dragged out of the 
disk (as seen in hydrodynamical simulations, e.g. Strick- 
land & Stevens (2000)). These dense clouds, although of 



low filling factor and again dominating the X-ray and opti- 
cal emission, would not necessarily give a limb-brightened 
optical or X-ray morphology. 

In principle, the filling factor of the X-ray emitting gas 
might be different in other parts of the flow. If density in- 
homogeneities in the flow are destroyed over a time scale 
less than or equal to the flow time, then the X-ray emis- 
sion in the outer parts of the wind might be dominated by 
a homogenized volume-filling hot component of the wind. 
Detecting filamentary or "clumpy" X-ray emission at lower 
surface brightness in the outer parts of the wind (with 
XMM-Newton or Chandra) would disprove this hypothe- 
sis. However, even if the X-ray emission appears relatively 
uniform the filling factor is not necessarily high. It is easy 
to conceive of situations where low filling factor clouds 
dominate the X-ray emission but whose spatial distribu- 
tion is so uniform as to appear smooth. 

It is only with the spatial resolution Chandra provides 
that it is possible to resolve (both spectrally and spatially) 
the inner outflow cones (within a few hundred parsecs of 
the starburst region) where unambiguous evidence of the 
origin of the soft X-ray emission in superwinds can be 
found. 

6. SUMMARY 

Chandra ACIS observations demonstrate that the soft 
X-ray emission in the less-obscured of NGC 253's super- 
wind's nuclear outflow cones is well coUimated, and ex- 
tremely similar in morphology and extent to the previously 
observed limb-brightened Ha outflow cone. This provides 
the flrst direct observational evidence that both X-ray and 
Ha emission come from low volume flUing factor gas, re- 
gions of interaction between a tenuous starburst-driven 
wind of SN-ejecta and the dense ambient ISM, and not 
from the wind itself. 

We measure the observationally and theoretically im- 
portant fllling factor of the X-ray emitting gas at between 
rj ^ A and 40 per cent. The observed morphology of the 
outflow also provides a purely geometrical upper limit of 
77 ~ 40 per cent. Only < 20 per cent of the observed 
X-ray emission can come from the volume-filling, metal- 
enriched, wind fluid itself. Spatially-resolved spectroscopy 
of the soft diffuse thermal X-ray emission reveals that the 
predominant source of spectral variation along the outflow 
cones is due to strong variation in the absorption on scales 
of ~ 60 pc, with little evidence for any change in the char- 
acteristic temperature of the emission. These observations 
are easily explained by, and fully consistent with, the stan- 
dard model of a superwind driven by a starburst of NGC 
253's observed properties. 

The important implication of this work is that we can 
not sec the majority the energy and gas (in particular the 
hot metal-enriched gas) transported out of NGC 253 by 
it's superwind, even in X-ray emission. If this is gener- 
ally true for all starbursts, the majority of IGM heating 
and enrichment by starbursts may be currently invisible in 
emission. This emphasizes the need for, and importance 
of, absorption line studies of warm and hot gas in these 
fascinating and important objects. 

We would like to thank the anonymous referee for their 
constructive criticism. The team that built and operates 
Chandra also deserves recognition, for their hard work. 
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and for making such an excellent telescope. This work was GO0-1008X. 
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